Ultrafast stimulated Raman parallel adiabatic passage by shaped pulses 
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We present a general and versatile technique of population transfer based on parallel adiabatic 
passage by femtosecond shaped pulses. Their amplitude and phase are specifically designed to op- 
timize the adiabatic passage corresponding to parallel eigenvalues at all times. We show that this 
technique allows the robust adiabatic population transfer in a Raman system with the total pulse 
area as low as 3 tt, corresponding to a fluence of one order of magnitude below the conventional 
stimulated Raman adiabatic passage process. This process of short duration, typically pico- and 
subpicosecond, is easily implementable with the modern pulse shaper technology and opens the pos- 
sibility of ultrafast robust population transfer with interesting applications in quantum information 
processing. 
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I. INTRODUCTION 

Complete or partial population transfer between 
ground states by pulsed laser fields is at the heart 
of many quantum processes [l| ranging from laser- 
controlled chemistry to modern quantum optics and 
quantum information, such as quantum gates [3|. 

Standard 7r-pulse techniques that require a specific 
pulse area are known to lack robustness 0, Adi- 
abatic passage, which allows the dynamics to follow a 
single eigenstate of the system, is often used to imple- 
ment such transfers due to its robustness with respect 
to field fluctuations and to the imperfect knowledge of 
the studied system. In particular, the stimulated Ra- 
man adiabatic passage (STIRAP) process, which induces 
a complete population transfer in a A system [6*1 , and var- 
ious extensions |7|, iSi, |9|] , have become very popular. The 
main advantage relies on its dark-state dynamics which 
does not involve the upper state in the adiabatic limit, 
and makes thus this process in principle immune to de- 
coherence. However this process requires in practice a 
large pulse area to satisfy the adiabaticity and prevents 
in most cases its use for pulses shorter than nanoseconds. 

On the basis of the adiabatic theorem, it is commonly 
believed that a large pulse area for adiabatic passage is 
inevitable and that it is the price to pay to get the ro- 
bustness. It has been however found that, for a two- level 
system driven by an appropriate shaped chirped pulse, 
the pulse area can be only typically I.Stt [l^, which is 
not much larger than the tt value (corresponding to the 
smallest value leading to a complete population transfer 
for a pulse of a laser perfectly tuned to the resonance. 
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the TT-pulse [Tl|). This finding has been made from a 
geometric picture showing the surfaces of eigenenergies 
as functions of the field parameters the paths that 
optimize the non-adiabatic losses are level-lines in the 
diagram of the distance between the two surfaces and 
correspond to instantaneous eigenvalues that are parallel 
at all times. 

In a A-system (where we assume that the two ground 
states are not coupled, e.g. for symmetry reasons), the 
naive simplest way to completely transfer the population 
between the two ground states is two successive tt pulses, 
one for each transition, leading to a total pulse area of 
2tt. One can also consider two exactly overlapping fields 
also giving a total pulse area of 27r [l| . It has been shown 
that the minimum area for such a process is VSn, which 
corresponds to the singular- Riemannian geodesic [ll|. 

We show in this paper that one can extend the opti- 
mization procedure making the dynamics follow parallel 
eigenvalues at all times to a A-system driven by two ap- 
propriately shaped fields of chirped frequencies. One gets 
a complete population transfer for a total pulse area as 
low as Stt with the standard robustness of adiabatic pro- 
cesses. This opens the possibility to implement ultrafast 
Raman adiabatic passage technique, i.e. in the pico-and 
subpicosecond regime, using the state-of-the-art technol- 
ogy of shaped femtosecond laser pulses (see for instance 
[l3l [T^ and its recent use and proposals for strategies of 
control [H, [H, [13, [S d ) . 

The paper is organized as follows. We first describe 
the technique of parallel adiabatic passage in A-systems. 
We next compare it with respect to the conventional 
STIRAP. Robustness with respect to fluctuations is ana- 
lyzed. Before concluding, we propose an implementation 
of the parallel STIRAP using pulse shaping techniques. 
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II. PARALLEL ADIABATIC PASSAGE IN 
A-SYSTEMS 

A. The system 

The A-system, with the two ground states denoted re- 
spectively |1) and 1 3), driven by two fields, respectively 
the pump and the Stokes fields, is described in the reso- 
nant approximation [l| in the dressed-state basis by the 
Hamiltonian 
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The available parameters that can a priori vary during 
time are flp and Qs (chosen real for simplicity) the pump 
and Stokes Rabi fi'equencies coupling the transitions |1) — 
1 2) and |3) — 12) respectively, A and d respectively the one- 
photon detuning (with respect to the pump) and the two- 
photon detuning: A = [^2 — Wi— Wp, 6 = UJs — LUi—LjJp+UJS 

with hxjj, j = 1,2,3, the energies of the corresponding 
state |j) and ujp, ujs the frequencies of the pump and 
Stokes fields respectively. 



B. Construction of parallel eigenvalues 

We apply the analysis that was made in Ref p^] for 
a two-state system to the three-state model: We require 
an adiabatic passage process such that the eigenvalues 
stay parallel at each time in order to optimize the adia- 
batic passage. We first remark that this imposes fields of 
chirped frequency. The use of (non-optimized) linearly 
chirped fields has been proposed in ^] and studied in 
[2l| with the aim of minimizing the intermediate-level 
population. We first determine sufficient conditions to 
achieve the parallelism of the eigenvalues and next give 
a time-dependent realization that can be implemented 
using modern pulse shapers. 

Adiabatic passage processes require the knowledge of 
the eigenelements. For the case of a two-photon detun- 
in g, th ey have been determined and studied for instance 
in |22l |. Here we require the parallelism of the eigenvalues. 
More precisely, denoting cj_,o,+ the three eigenvalues sat- 
isfying at each time uj^ < ujq < uj+ and \ip-,o,+) the 
corresponding eigenstates, we impose for reasons of sym- 
metry a;+ — luq = LUo — lu- and that \ipo) transports adia- 
batically the population, thus connecting at early times 
the initially populated state |1) and at late times the tar- 

t — ^-f oo 



get state |3): |1) <— 
are thus of the form 



|3). The eigenvalues 
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and the condition 



= [9{nl + nl)~4:(2A~S){2S- A) 
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The initial connection, when Qp = Cls ~ 0, imposes 
for the time-dependent detuning 6(~oo) ~ — A(— oo) = 
rio/2 corresponding to uji^i—Qo) — 0, and cj±(— oo) — 
±r2o/2. The final connection can be chosen in two ways: 
Either ti;_|_(-|-oo) = 0, ti;o(-|-oo) = — Slo/2, a;_(+oo) = 
— ilo, or ci;+(-|-oo) = fio, a;o(+oo) = cl'_(-|-oo) — 0. 

We have chosen below the latter situation correspond- 
ing to i5(-l-oo) = Oo/2 and A(-l-oo) = f^o. At this stage, 
we expect many solutions for the four available param- 
eters from Eqs. ^ and (|3|) with boundary conditions: 
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C. Numerical implementation 
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FIG. 1: Transition Probability for a level line dynamics as 
function of time for Qq = 5.8/T (corresponding to the to- 
tal pulse area A « S.Bvr). Upper frame: Population Pj = 
|(j|(/)(t))P with (j){t) the state solution. Middle frame: The 
field parameters (in units of 1/T) A{t), S(t) given by Eqs. ([5} 
(for Of = 0), and ^lp{t) and ils{t) determined from Eqs. (|3]) 
and ((3}. Lower frame: The eigenvalues (in units of 1/T). 



To find a simple implementable solution, we addition- 
ally impose the time parametrization as a monotonic 
smooth function for A(i), symmetric around t = 0, and 
a symmetric Gaussian shape for 6{t): 
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with a and f3 to be chosen. The quantity T corresponds 
to a characteristic time of the duration of the process. 
The choice a = is simple since it leads to a constant 
two-photon detuning; it will be shown to give already 
very good results. For given a and /3, the choice fio comes 
down to choose a unique level line for the problem and 
set the values of the Rabi frequencies rip(t) and ris(t) 
from Eqs. (3) and (4). Note that the parametrization 
([5]) is somewhat arbitrary. Other smooth shapes for the 
detunings lead to very similar (but not strictly identical) 
results. 

We define the total pulse area for the process as 



A= dt^n%{t) + nl{t). (6) 

J —oo 

This characteristic quantity is known to be at least tt 
when one considers a two-state transition (for instance 
the 1-2 transition by the pump field) The pulse 

area has a clear-cut dynamical meaning only for resonant 
pulses in two-level systems. In other (nonresonant, multi- 
level) situations, different pulses with the same individual 
area can produce different results even in near-adiabatic 
cases. This choice of the total pulse area takes also into 
account the delay between the pulses. Considering for 
simplicity two identical delayed pulses of individual area 
Ap, the area A is bounded between \/2Ap (no delay) 
and 2Ap (consecutive pulses without overlapping). This 
criterion favors thus overlapping pulses since they lead to 
faster processes. In the A system, the use of two succes- 
sive TT pulses, respectively for the pump and Stokes fields, 
or of two overlapping fields corresponds both to A — 2tt 
[l| . However it is known that the minimum area for such 
a process is VStt [Ti[. An alternative quantity of inter- 
est is the total fluence corresponding to an integrated 
intensity: 

/ + 00 
dt{n%{t) + nl{t)). (7) 
-oo 

Figure [1] shows the resulting dynamics corresponding 
to parallel eigenvalues with fto — 5.8/T and a = 0. We 
obtain A ~ S.Gtt and « AO/T in these conditions. 
The obtained final population transfer P3(-|-oo) « 0.995 
shows a very efficient population transfer by adiabatic 
passage for a quite modest pulse area, only approximately 
twice the minimum required area y/Sn. One can no- 
tice the counterintuitive sequence of the pulses (i.e. the 
Stokes field coming first before the pump field), which 
appears here as a consequence of the level line procedure 
for the choice ([5]). We notice a non-negligible transient 
population in the upper state 2 since here the dynamics 
is not governed by a dark state. We however anticipate 
that a fast enough process, typically subpicosecond, will 
allow one to neglect the resulting loss from the upper 
state. 



III. COMPARISON WITH THE 
CONVENTIONAL STIRAP 
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FIG. 2: Upper frame: Final population transfer for the par- 
allel STIRAP with a = (full line), linearization of A{t) [see 
Eq. ITO])] (large-dotted line) , the conventional STIRAP with 
the delay r — I.IT (dashed line) and with the delay t = T 
(dot-dashed line), and parallel STIRAP with a = O.l/T and 
l3 — 1.25 (thin-dotted line) as a function of the total pulse 
area A. Middle frame: logarithmic scale of the deviation 
from the complete population transfer. Lower frame: The 
corresponding QqT for the parallel STIRAP's. 
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FIG. 3: Same as Fig. [2] but as a function of the fluence (in 
units of 1/T) in logarithmic scale. 
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Figures [2] and [3] show a comparison of the Stimulated 
Raman paraUel adiabatic passage using ([5]), for a con- 
stant {a = 0) and a non-constant (a 0) two-photon 
detuning, with the conventional STIRAP in terms of 
population transfer efficiency with respect to, respec- 
tively, the pulse area A and the fluence J-. The con- 
ventional STIRAP process is considered with a coun- 
terintuitive sequence of Gaussian pulses of same peak 
values that are delayed of t: ftp — ilmaxe"-*^^^^''^"^' , 
= f^maxe"'^'"''^''^^/"^' . Two values for t are consid- 
ered (with a difference of 10%) in Figs. [2] and H We 
can first notice a significant dependence on this delay for 
the relatively small areas that are here considered and 
thus a relative lack of robustness. Note that we have also 
studied the influence of the delay between the two pulses 
for the parallel STIRAP: We have obtained a much less 
significant dependence. 

For the constant two-photon detuning case {a — 0), 
the population transfer to state 3 by the parallel strat- 
egy is already very efficient for A~ 3.6n (corresponding 
to the dynamics considered in Fig. [T]) and gets monoton- 
ically better for higher values of A. This value A « B.Gtt 
corresponds to a weak efficiency for the conventional STI- 
RAP. We recover this efficiency for the conventional STI- 
RAP for a pulse area approximately three times larger 
than for the parallel strategy, corresponding to a fluence 
approximately one order of magnitude larger (see Fig. 
[3]). This shows the significant superiority of the paral- 
lel strategy with respect to the conventional STIRAP. 
The deviation from the complete population transfer dis- 
played in Figs. [2] and [3] with a logarithmic scale shows 
the high efficiency of the parallel strategy. We obtain a 
transfer with an accuracy to more than 4 digits for an 
area larger or equal to only 5.57r. We remark that such 
accuracy for the population transfer is not possible for 
the conventional STIRAP (see also [23]). 

Figures [2] and [3] shows that the efficiency of the paral- 
lel strategy can be improved using a non-constant two- 
photon detuning (a ^ 0): We obtain the same transfer 
(considering a deviation from the complete population 
transfer not smaller than 10~^) with a significant smaller 
area for the non-constant two-photon detuning case with 
respect to the constant case. Numerics shows that the 
transfer is generally better when a and /3 are such that 
the two pulses approximately overlap (requiring a 0). 
The transfer efficiency is shown in Fig. [^in that case for 
a = 0.1 ilo and /3 — 1.25. We obtain a very good transfer 
efficiency already from A ~ 37r. The gain for the fluence 
with respect to the conventional STIRAP is more than 
one order of magnitude as shown in Fig. [31 



mainly based on the use of the Davis-Dykhne-Pechukas 
formula in the time complex plane with analytic func- 
tions as the pulse parameters |10l |. 

We consider a white noise of width F = G.5T on the de- 
tunings Apit) and Spit), on the Rabi frequencies flp_j{t), 
j = P, S, and on their areas, defining 
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with ilj (t) the Rabi frequency of the field j without fluc- 
tuation, < Aj{t) < 1 the shape of the fleld j without 
fluctuation, A(i) and S{t) the instantaneous detunings 
without fluctuation, and —0.5 < Vk < 0.5, k = 1, • • • ,4, 
four (independent) random numbers. The random num- 
ber ri is drawn once for each realization of a complete 
dynamics. It allows one to vary randomly the pulse ar- 
eas. The random numbers r^, k > 2 have been drawn at 
each time step of calculation (here taken as r/300) for 
all realizations. This allows one to model instantaneous 
fluctuations of the detunings and of the pulse shapes. We 
have determined the instantaneous populations shown in 
Fig. |4]by averaging over many realizations of time histo- 
ries (Monte Carlo simulation) . The averaged final popu- 
lation transfer is here P3{+oo) w 0.993 (to be compared 
to P:j{+oo) Ri 0.996 without fluctuations). 
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IV. ROBUSTNESS OF PARALLEL STIRAP 

The robustness of the process with respect to an im- 
perfect knowledge of the pulse areas is shown in Fig. [21 
Sources of other fluctuations can be various and the ro- 
bustness of the parallel strategy is questionable since it is 



FIG. 4; Same as Fig. [Tj but for flo = 5.4/r, a = 0.1, and 
l3 = 1.25, and with fluctuations for the amplitudes and instan- 
taneous frequencies. The field parameters and eigenvalues are 
shown for a single realization. 



V. IMPLEMENTATION OF PARALLEL STIRAP 
BY PULSE SHAPING TECHNIQUES 

The implementation of the paraUel STIRAP requires 
the design of the pulses and frequencies using the state- 
of-the-art technology of pulse shaping, such as for in- 
stance liquid cristal spatial light modulators. Two such 
modulators allow the independent shaping of the spec- 
tral phase and amplitude of an input femtosecond laser 
pulse [l^, [13] ■ We consider here for simplicity a constant 
two-photon detuning (a = 0). In that case, we have to 
chirp both field frequencies exactly in the same way. The 
difference between the two fields is their mean frequency 
and their delayed envelope. We can determine approxi- 
mate curve fitting of the envelopes for a = 0: 



\ (9) 



with -I- (— ) holding respectively for the Stokes (pump) 
field. Figure [5] shows, for these fields of phase correspond- 
ing to (O with a = (cf. conditions of Fig. 1), their cor- 
responding spectral amplitudes and phases, respectively 
denoted A[uj) and and defined through the complex 
field as £'A(t)e*'^(*) = ^ / dw^(w)e^['^*+'*(")l (for each 
field). Here (p{t) is the instantaneous phase, leading the 
instantaneous frequency dip/dt, < K{t) < 1 the shape of 
the field amplitude, and £ the peak amplitude of the con- 
sidered field. The instantaneous frequencies of the pump 
field LOp = dipp/dt and of the Stokes field lus = dips/dt 
are connected to the one- and two-photon detunings as 
follows: u!p(t) = uj2—uJi—A(t), ujsit) = uj2—uj3+S—A{t). 
We have chosen for convenience the mean frequency of 
the pump as cjq.p = 1^2 — ti^i — r2o/4 = ^p(O) and of 
the Stokes as wq.s = ^2 — 1-^3 + S — rio/4 = ws(0). We 
have considered fields of full width at half maximum 
TpwHM — 500 fs corresponding to Oq = 17.1 THz for 
the considered case Hq = 5.8 /T. For typical transition 
strengths of 1 Debye, this leads to fields of peak intensi- 
ties 15 GW/cm^. For each field, the spectral parameters 
of Fig. [5] have been generated from a Gaussian Fourier 
transform limited pulse centered at the mean frequency 
luq with a full width at half maximum of 100 fs, shaped in 
phase and amplitude using two spatial light modulators. 
Note that we have considered the realistic effect of pix- 
elization for the modulators (using 320 pixels, almost not 
noticeable at the scale of Fig. O which does not show 
any noticeable difference for the output pulses and the 
corresponding dynamics of Fig. [TJ 



For practical implementation, it could be of interest to 
alternatively use, with the fitting Gaussian pulses ([HI), lin- 
early chirped frequencies constructed from the lineariza- 
tion of ([5]) around < = 0: 
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FIG. 5: Tailored spectral normalized amplitudes A (upper 
frame) and phases (lower frame) for the pump (full lines) and 
Stokes (dashed lines) fields from 100 fs Fourier transform lim- 
ited pulses (upper frame, dotted line) of mean frequency cjq 
to produce the fields of Fig. [1] 



is maximum. This would indeed allow 
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(10) 



where ^/rip 

one to use for instance gratings to produce the chirp in 
stead of modulators (see for instance [2j] ) . In this case, 
the duration of the pulse is connected to the slope of the 
chirp and the pulses have to be delayed additionally. Its 
efficiency is shown in Figs. [5] and [3J It is not as effi- 
cient as the level line strategy but is still well superior to 
the standard STIRAP. For instance, one can see that the 
deviation from the complete transfer is less than 10~^ 
with a pulse area larger than S.Stt for this linearization, 
corresponding to a fluence larger than lO^'^/T. One re- 
covers this accuracy for STIRAP for a fluence larger than 
IO.Stt, corresponding to a fiuence larger than lO^ '^/T (al- 
most one order of magnitude larger). This linearization 
also allows an error smaller than 10"'' for an area larger 
than IOtt (it should be larger than only S.Stt for paral- 
lel STIRAP). Such accuracy is not reachable using the 
conventional STIRAP. 



VI. CONCLUSION AND DISCUSSIONS 

In conclusion, we have developed a strategy of adia- 
batic passage in A-systems with field designed such that 
the eigenvalues stay parallel at each time. This strat- 
egy allows one to combine the energetically efficiency of 
"TT-pulse" and related strategies with the robustness of 
adiabatic passage. It has been shown in particular to 
be superior to the conventional STIRAP. It is expected 
to be easily implementable using modern tools of pulse 
shaping. 

To implement the proposed strategy in real systems, 
one has to consider several limitations due to the pres- 
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ence of other states and to lossy processes such as ioniza- 
tion. In order to avoid the latter, the requirements are 
to consider a system which does not show a one-photon 
resonance into the continuum from the excited state and 
to use moderate field intensities, typically not more than 
100 GW/cm^ for atoms. Additional non-resonant states 
lead to Stark shifts which can be considered in general 
negligible with respect to the one-photon resonances con- 
sidered here. However for the considered field intensities, 
they can play a role in deforming the dressed energies. In 
that case one has to recalculate the field parameters such 
that the eigenenergies stay parallel. The conditions that 
permit these additional states to stay non resonant can 
be estimated as A^-'^T > 1 and < A\^^ with A[:'^ 
the one-photon detuning for the state j and the field k 
and fl^^^ the corresponding peak Rabi frequency. 

This parallel adiabatic passage strategy should be 
adapted to produce in an optimal way superpositions of 



states such as in the case of fractional STIRAP [Ssj . This 
should find applications in quantum information process- 
ing, for instance to implement fast quantum gates or 
quantum algorithm (such as the implementation of the 
parallel adiabatic passage for the quantum search [2^). 
The high efficiency of the parallel strategy shown in Fig. 
2 is compatible with such applications imposing severe 
limitations on the admissible error, typically on a level 
below 10"*^ (see, e.g., [1^ and references therein). 
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